Transport properties of a superconducting single-electron transistor 
coupled to a nanomechanical oscillator 
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We investigate a superconducting single-electron transistor capacitively coupled to a nanome- 
chanical oscillator and focus on the double Josephson quasiparticle resonance. The existence of two 
coherent Cooper pair tunneling events is shown to lead to pronounced backaction effects. Measuring 
the current and the shot noise provides a direct way of gaining information on the state of the oscil- 
lator. In addition to an analytical discussion of the linear-response regime, we discuss and compare 
results of higher-order approximation schemes and a fully numerical solution. We find that cooling 
of the mechanical resonator is possible, and that there are driven and bistable oscillator states at 
low couplings. Finally, we also discuss the frequency dependence of the charge noise and the current 
noise of the superconducting single electron transistor. 

PACS numbers: 85.85.-|-j,73.23.Hk,73.50.Td 



I. INTRODUCTION 

The cooling of nanomechanical systems by measure- 
ment has received a lot of attention recently. Various 
procedures like the laser sideband cooling schemes 
developed for trapped ions and atoms fi. have been 
proposed as ways to significantly cool a nanomechanical 
resonator (NR) coupled to a Cooper-pair box,— i^i^i^ a flux 

qubitf^ii quantum dotsjii, trapped ions^ and optical 
cavities ,i3-i4. is. 16. 17. is. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34 

On the experimental side, optomechani- 
cal cooling schemes have been shown to be 
promisingii^iiiii^ii^iilii^ii^i^a the NR was cooled to 
ultra-low temperatures via either photothermal forces or 
radiation pressure by coupling it to a driven cavity. 

Another important nanoelectromechanical measure- 
ment device which both holds the possibility of very 
accurate position measurements^^ as well as of cooling 
of an NR, is a superconducting single-electron transis- 
tor (SSET). Shortly after the theoretical proposals pre- 
dicting the potential of the SSET to cool a nanome- 
chanical systemjSiiSiiS this effect has been experimen- 
tally observediiS Using other detectors for NRs such 
as normal-state single-electron transistors^^ or tunnel 



junctions 



37.68 



it is very difficult to cool the nanomechan- 



ical system or drive it into a non-classical state. These 
detectors usually act as heat baths with effective tem- 
peratures proportional to the transport voltage, which 
is in practice higher than the bath temperature. The 
SSET system, on the other hand, shows sharp transport 
resonances. At those the effective temperature is voltage- 
independent and can be made very low. 

To achieve such challenging goals as ground-state cool- 
ing of NRsfi^ or the creation of squeezed oscillator states, 
a better understanding of the transport properties of the 
coupled SSET-NR system is required. This system is 
schematically shown in Fig.[T] Depending on external pa- 



rameters such as the gate voltage Vq, the bias voltage V, 
and also the superconducting gap A, the SSET supports 
different types of resonance conditions. The two most 
prominent ones are the so-called Josephson quasiparticle 
(JQP) and the double Josephson quasiparticle (DJQP) 
cyclej^i^S Whereas the former involves the coherent tun- 
neling of a Cooper pair at one of the two junctions fol- 
lowed by a successive tunneling of two quasi-particles at 
the other junction, the latter consists of four steps (illus- 
trated in Fig.[2]below) that involve a Cooper pair tunnel- 
ing at each of the junctions and a quasi-particle tunnel- 
ing at each of the junctions. The transport properties of 
the SSET coupled to an NR close to the JQP resonance 
have been analyzed in a recent theoretical work4S Here, 
we focus on the analysis of the same coupled quantum 
system at the DJQP resonance. Since the JQP is a one- 
dimensional resonance in the parameter space spanned by 
Vg and V and the DJQP is a zero-dimensional resonance 
in the same parameter space, all action and backaction 
effects close to the DJQP resonance are much more pro- 
nounced than close to the JQP resonance. This is of 
crucial importance if one wants to manipulate the state 
of the NR by measurement of the SSET detector be- 
cause, in experiments, the typical coupling between the 
two quantum systems turns out to be rather weak. 

We analyze how the NR can be cooled below the tem- 
perature of the external heat bath and how it can be 
brought into a (non-thermal) driven state at the DJQP 
resonance. Under certain conditions, we find signatures 
of bistable solutions of the coupled quantum system of 
NR and SSET. It is of particular interest and experi- 
mental relevance, to know how a successful cooling of the 
NR or the preparation of a driven state can be observed 
in transport properties of the SSET such as its current 
or current noise. We show that there is a one-to-one 
correspondence between interesting state preparations of 
the NR and the transport properties of the SSET. This 
provides a powerful and feasible tool to initialize and ma- 
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FIG. 1: Schematic setup of the SSET- resonator system; Two 
superconducting leads at vohages Vl and Vr are coupled by 
tunnel junctions to a superconducting island. Its chemical po- 
tential can be tuned by a gate voltage Va- A nearby nanome- 
chanical oscillator acts as an s-dependent gate. 



nipulatc NR quantum states by measurement. 

The article is organized as follows. In section |lTl 
we present the model for the coupled quantum system 
of NR and SSET, discuss the different approximation 
schemes of the analytical solutions as well as the cal- 
culation scheme behind the exact numerical solution of 
the underlying master equation. Then, in section ITTTl we 
analyze the oscillator properties by means of the differ- 
ent methods, identifying interesting quantum states of 
the NR due to its coupling to the SSET. Subsequently, 
in section IIVI we discuss the current of the SSET de- 
tector and, in section |Vl the charge and current noise. 
It turns out that the combination of the two transport 
properties is sufficient to clearly identify a successful cool- 
ing or driven-state preparation of the oscillator. Finally, 
we present our conclusions in section IVII Details of the 
calculations arc contained in the Appendices. 



II. MODEL 

The system under investigation consists of a supercon- 
ducting single-electron transistor (SSET) which is capac- 
itively coupled to a nanomechanical resonator (NR) as 
shown schematically in Fig. [T] The total Hamiltonian of 
the system reads 



H = Hr 



Hn 



Hj + Ht + He + 



N 



N,I 



(1) 



The first three terms Hl h j are standard BCS Hamil- 
tonians and describe two superconducting leads (left and 
right) and a superconducting island, 



Hn 



(2) 



Here, Caka arc annihilation operators for quasiparticles of 
momentum k and spin a in the system a {a = L,R,I). 
The dispersion relation Caka- accounts for the supercon- 
ducting gap of width 2A which we assume to be equal 
for the three systems. The chemical potentials in the left 
and right leads are determined by the applied bias volt- 
age V = Vl — Vr, while the island chemical potential can 
be tuned by applying a gate voltage Va (see Fig. [IJ. 



The left and right leads are connected to the central 
island by quasiparticle tunneling and Cooper pair tun- 
neling. Denoting by (f>a the superconducting phase dif- 
ference at the junction a = L, R, we use the following 
quasiparticle tunneling term 



H' 



T,qp 



E 



-L.R 



E 

k.q.a 



(3) 



where T^q are the tunneling amplitudes which can be 
used to calculate^ the quasiparticle tunneling rates Fi^i?. 
Cooper pair tunneling is accounted for by the term 



T,CP — 



E 

i=L,R 



J a cos ( 



(4) 



where Jq are the Josephson energies of the two junctions. 
Hence, the total tunneling Hamiltonian is given by Ht = 

Hr.qp + Ht,cp- 

The final ingredient for the SSET Hamiltonian is the 
Coulomb energy of the island. If we denote by and 
nfl the number of electrons that have tunneled from the 
island to the left and right lead, respectively, then n = 
—ul — nn is the excess number of electrons on the island. 
The charging term can be written as 



Hc^Ecin + nof + eVuR, 



(5) 



where Ec is the charging energy and no can be controlled 
by the gate voltage (see Appendix In terms of the 
capacitances of the two junctions Cl,r, the gate Co and 
the resonator Cn, the charging energy is given by Ec = 
e2/(2Cs), where Cs = Cl + Cr + Cg + Cn is the total 
capacitance. 

Next, we focus on the coupling of the SSET to the 
NR. The latter can be regarded as a harmonic oscillator 
of frequency fl and mass M and is therefore described by 
the Hamiltonian 



H, 



N 



r 

2M 



(6) 



The NR is held on a constant voltage Vm and hence acts 
on the SSET as an additional gate with an a;-dependent 
capacitance Cn{x). Therefore, the presence of the NR 
modifies the charging term He- Expanding the contri- 
bution for small displacements x and retaining only the 
lowest order, one finds that the coupling between SSET 
and NR is given by 



H. 



N.I 



= -^5 



(7) 



where the coupling constant A depends in a non-trivial 
way on the voltages and capacitances of the system 
and can be regarded as an effective parameter (see Ap- 
pendix]^. Note that this expansion is only valid for dis- 
placements X which are small compared to the distance d 
between the SSET and the NR, i.e. x/d 1. Upon con- 
tinuing the expansion, one encounters terms proportional 
to and to which will be neglected here. 
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FIG. 2: Illustration of the DJQP cycle; (i) Cooper pair 
tunneling through the left junction, (ii) quasiparticle tunnel- 
ing through the right junction, (iii) Cooper pair tunneling 
through the right junction and (iv) quasiparticle tunneling 
through the left junction. 



Due to the complex structure of the full Hamilto- 
nian ([T|) one should not hope for an exact solution in 
all regimes. Instead, we will make several assumptions 
which will enable us to investigate the transport proper- 
ties of this system at a particular point in the parameter 
space. 

First, wc will briefly review the transport properties of 
the bare SSET without coupHng to the NR. While the 
capacitances, Josephson energies and quasiparticle tun- 
neling rates are essentially determined by the experimen- 
tal setup, the most important tunable parameters are the 
bias voltage V and the gate voltage Vg- The transport 
properties of the SSET are then determined by how these 
voltages are related to the superconducting gap 2A and 
the charging energy Ec- 

For high bias voltages eV > 4 A, the difference in chem- 
ical potentials allows quasiparticlcs on both junctions to 
overcome the superconducting gap and a quasiparticle 
current can flow. But even for lower bias voltages, one 
observes a finite current at certain values of the gate volt- 
age. A possible mechanism is the Josephson-quasiparticle 
(JQP) resonance which is a cyclic process that starts with 
the tunneling of a Cooper pair on one of the junctions fol- 
lowed by two subsequent quasiparticle tunneling events 
on the other junctioni^i^ This process is possible above 
a lower bias voltage threshold, eV > 2 A -|- Ec- 

For even lower bias voltages, isolated current reso- 
nances can be observed which are due to the onset of 
the double Josephson quasiparticle (DJQP) resonance. 
A schematic picture of this process is shown in Fig. [2] It 
starts with a Cooper pair tunneling across, say, the left 
junction. Next, a quasiparticle tunnels out through the 
right junction, followed by a Cooper pair. Finally, after a 
quasiparticle tunnels through the left junction, the initial 
system state is reached again. 

This process is energetically allowed only in a restricted 
parameter regime: Cooper pair tunneling is only possi- 
ble if the chemical potentials of the lead and the island 
(taking into account the Coulomb energy) are on reso- 
nance while quasiparticle tunneling requires a difference 



in chemical potentials sufficient to overcome the super- 
conducting gap. For the DJQP process, it is easy to show 
that the resonances occur at bias voltages eV = 2Ec and 
half-integer island charges jiq. 

The parameter regime which we investigate is there- 
fore characterized by a charging energy Ec, a supercon- 
ducting gap 2A and a bias voltage V which are of the 
same order of magnitude. Roughly speaking, these en- 
ergy scales are very large compared to the quasiparticle 
tunneling rates T^n, the Josephson energies Jl.r and 
the oscillator energy Q. 



A. Derivation of a Liouville equation 

Due to the small tunneling rates, only sequential tun- 
neling will contribute to the transport whereas higher- 
order (cotunneling) processes are suppressed. This sug- 
gests describing the system by a master equation in the 
Born-Markov approximation. 

For this purpose, we treat the BCS Hamiltonians 
Hl + Hn + Hi as a fermionic bath for the remaining 
system. Then, system and bath are only coupled by 
the quasiparticle Hamiltonian HT,qp- Using the Born 
approximation corresponds to disregarding cotunneling 
processes while the Markov approximation is valid as long 
as there is a separation of time scales between the system 
and the bath degrees of freedom. Introducing the system 
and bath Hamiltonians 



Hs — He + Ht.cp + H]^ 
Hb = Hl + Hji + Hi , 



H 



NJ 



(8) 
(9) 



and using the Born-Markov approximation leads to the 
following master equation for the reduced density matrix 
p{t) of the system. 



p{t) = Cp{t) 



(10) 



= --[Hs,pit)] 
1 

- J dTTTB[HT.,qp,[HT.,qpi-T),pit)(g, Pb]] , 

where pB is the bath density matrix. The time depen- 
dence of the Hx^qp operator is governed by the Hamil- 
tonian Hs + Hb- The density matrix p contains infor- 
mation only about the charge and the oscillator degrees 
of freedom and can, for example, be written in the ba- 
sis \n,nii,x) of island charge states \n), the amount of 
charge jrii^) which has tunneled through the right junc- 
tion, and the oscillator coordinate \x). This approach 
allows the calculation of the transport properties of the 
system via charge counting.— In order to investigate 
the transport at the DJQP resonance, it is sufficient to 
consider a finite number of basis states for the island 
charge n. As a single DJQP cycle involves four charge 
states, we can restrict the basis to the states | — 1), |0), 
1), and 1 2) which significantly reduces the complexity of 
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the problem, since it is thus sufficient to study a reduced 
density matrix as described in section III CI This choice 
of charge states corresponds to no = —1/2. 

As long as one is only interested in oscillator properties 
or the current through the SSET, the states can be 
traced out, and an effective master equation acting on 
the Hilbert space of island charge and oscillator posi- 
tion, spanned by the states \n,x), can be obtained. On 
the other hand, for the calculation of the current noise, 
the jrifl) degree of freedom has to be taken into account 
explicitly, as will be explained later on. For now, we 
proceed with the rijj-independent case. 

As in the case of the JQP^ the Liouvillian obtained 
from Eq. (fTO|) can be written as a sum of three contribu- 
tions 

^ = ^Hs + ^qp + ^CL , (11) 

where Chs governs the coherent evolution of the system, 
Cqp is a dissipative term due to quasiparticle tunneling 
and CcL is a Caldeira-Leggett type contribution intro- 
duced to model the coupling of the harmonic oscillator 
to a finite-temperature environment. Explicitly, 

ChsP^-'j:[Hs,p], (12) 
^qpP = ^LP~iflPP^-ifl - _i, p} (13) 

C-cLP^ -^[xA^.p]] - ^^''^^\ x,{v,p)] , (14) 

where p^j = \j) {k\ acts on the charge states of the island 
and is utilized here to describe the quasiparticle tunneling 
event that changes the charge state from \k) to The 
energy dependence of the quasiparticle tunneling rates 
Fl and F^ is weak and will therefore be neglected in 
the following. The diffusion constant D and external 
damping rate 7cxt are related via a fluctuation-dissipation 
relation 

i^^M,,.,-coth(^^j 

= Af7cxt hn (^{nosc) + ■ (15) 

For hfl <C kBTs, the diffusion constant can be approxi- 
mated by £) = M-fcxtkBTs- 

In the general case, it is neither possible to calculate ex- 
actly the steady-state properties nor the transport prop- 
erties of the coupled system using the Liouville super- 
operator (jlip . Thus, approximation schemes must be 
employed. In the following two subsections, we describe 
in detail the two complementary approximation schemes 
we use to study the coupled SSET-oscillator system. 



B. Mean- field approach 

Physical quantities can be calculated by evaluating the 
matrix elements of the density matrix p{t). It turns out 
that oscillator properties and the average current can be 
written in terms of expectation values of the form 

{x^v"'pkj) = Tiw Yl "^1 2;"«"V(0 \J, riR + k^ j) , 

(16) 

where x and v are the position and velocity operators of 
the oscillator and pkj = \j) (k\. The trace over the oscil- 
lator degrees of freedom Trosc will be used in the position 
basis Trosc (•) = / dx {x\-\x) as well as in the phonon num- 
ber basis where Trosc (•) = I]^,,=o ("osc| ■ |"-osc). For the 
uncoupled SSET tuned closely to the DJQP resonance, 
the average current can be calculated straightforwardly, 
as taking the matrix elements of the master equation (|10p 
leads to a closed set of equations4^ If the NR is included, 
however, the coupling terms will lead to equations in- 
volving matrix elements of the form (xpkj ) . Calculating 
their time evolution leads to ever higher-order terms of 
the form {x"v"^pkj) , so that the set of differential equa- 
tions never closes. Hence, a truncation scheme is needed. 
A standard route is to truncate the system of equations 
by assuming a vanishing nth-order cumulant ((x'^pkj))- 
This allows one to rewrite nth-order expectation values 
in terms of expectation values of order n — 1 and hence 
to arrive at a closed, albeit non-linear, set of equations. 

We use and compare these approximations for n ~ 
1 (which we call thermal-oscillator approximation) and 
n = 2 (Gaussian approximation). These two levels of 
approximation are related to what was called 'mean 1' 
and 'mean 2' in Ref. In order to give an estimate 

of the physical quality of the truncation scheme, we also 
compare our results to exact numerical calculations. 

Whereas the expectation values of the form are 
sufhcient for the calculation of the oscillator properties 
and the SSET current, the calculation of the noise re- 
quires a slightly extended approach. In order to keep 
track of the transfered charge, one has to investigate the 
dynamics of n ^-resolved expectation values. It turns out 
(see Appendix[E| that the noise can be rewritten in terms 
of expectation values of the operators 

Ptf = \J) {k\ ® Wh) KI . (17) 

Note that only elements which conserve the number of 
charges, j + ^ k + n'j^, are finite. An analogous 
truncation scheme can be applied to expectation val- 
ues containing these operators. Similar approaches have 
been used extensively to describe nanoelectromechanical 
systems^^^^^^ 

C. Numerical solution of the Liouville equation 

To complement the analytical mean-field approach de- 
scribed in the last subsection, we also use a numerical 
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approach to study the properties of the NR coupled to 
an SSET near the DJQP resonance. First, we present 
the approach taken for the calculation of the current and 
the oscillator properties, where the nn degree of freedom 
plays no role. Subsequently, we will demonstrate how to 
extend this approach for the noise calculation, where tir 
has to be taken into account. 

To calculate the current and the oscillator properties, 
we write the density matrix in the |n, riosc) basis, with 
JT-osc being the phonon quantum number of the oscilla- 
tor and n the charge of the SSET. The spectrum of the 
harmonic oscillator is naturally not finite, so we need 
to truncate it and consider only its N„iax lowest energy 
eigenstates. To describe the DJQP cycle, the reduced 
density matrix p is of dimension {4:Nmax x ^Nmax)- 

The Liouvillc supcropcrator £ is a two-sided operator, 
in the sense that it acts both from the left and the right 
of the density matrix [cf. Eq. (fT2|) ] . It can be transformed 
to a single-sided operator using a property of the matrix 
vectorization operation: the vectorized form of a product 
of three {4:Nmax x 'iNmax) matrices A,B,C can be written 
as a single product of an {16N^^^ x 16N^„^^) matrix with 
an (16iV,^^^ X 1) vector via the relation vec(ABC) = 
(C"^ (8> A)vec(B), where the superscript T denotes the 
matrix transposition and ® a Kronecker product.— The 
matrix representation of the Liouville superoperator is 
therefore of order (16iV^„^ x 16N'^^^^). 

To illustrate how the aforementioned vector identity 
can be used, we apply it to the coherent evolution con- 
tribution to the Liouville equation [Eq. p^ ]. In this 
case, we find 

vec{CHsP) = -^{l® Hs + ®1) vec(/9), (18) 

where the matrix representation of the identity matrix I 
and of the system Hamiltonian Hs is of order {4:N„iax x 

To find the vectorized form of the stationary density 
matrix pstat, defined from Cpstat = 0, we calculate the 
null-space of the Liouvillc matrix. Using the normaliza- 
tion condition Tr[/9stat] = 1, the stationary density ma- 
trix Pstat can be determined uniquely. The bad scaling 
^{■^max) the Liouvillian size with the truncation point 
in the oscillator spectrum makes the numerical eigenvalue 
problem very challenging. Luckily, in this problem the Li- 
ouville matrix displays a high sparsity degree, and sparse 
eigensolvers can be used. Our implementation uses the 
shift-invert mode of the ARPACK^ eigensolver in com- 
bination with the PARDISO^i^ linear solver to com- 
pute the first few {^^ 5) eigenvalues of C with the lowest 
magnitude (Aq) as well as the associated eigenvectors. 
The calculated magnitude of the smallest eigenvalue can 
be used to verify the validity of the truncation scheme: 
when enough Pock states are kept we find |Ao| < 10~^'^ 
which is below the desired precision limit of 10~^°. 

To improve the speed of the calculation and, more 



importantly, to increase numerical accuracy, we do not 
need to explicitly solve for those matrix elements of Pstat 
which, due to the considered Hamiltonian, have to be 
zero. For example, coherence can only be created be- 
tween two charge states |fc) and |j) if |fc — j| = 2, 
since only these pairs of states arc coupled by Joseph- 
son tunneling. Therefore, all density-matrix elements 
{k\p\j) where \k — j\ is odd are zero. Using this ar- 
gument, the size of the Liouvillc matrix can be reduced 

to mlax X SA^rLJ- 

The use of sparse solvers also minimizes the required 
memory for the calculation of the eigenvalues, allowing 
problems of relatively large size {N^ax ^ 150) to be 
solved on a desktop computer. Also, we note that, con- 
trary to what was discussed in Ref. [S^] , no manual pre- 
conditioning was needed to achieve high numerical accu- 
racy. To allow for the numerical approach to be used in 
the driving regime, where we expect the average energy 
of the oscillator to be relatively high, we had to make a 
supplementary approximation. In this case, we assumed 
that coherence could develop only between states of the 
oscillator that arc not too far away in energy from each 
other, setting (nosc|p|"o.sc) = for l"osc - «oscl ^ 60. 
This allowed for N^ax to be set as high as 750 on a stan- 
dard workstation. Moreover, in the cases where it was 
possible to compare directly the results of the calcula- 
tion with and without this last approximation, we found 
that they were identical (within our numerical accuracy). 

While this approach is viable for the calculation of the 
oscillator properties and the current, it fails to keep track 
of the tunneled charge n^j and thus cannot be used to 
calculate the current noise. A straightforward inclusion 
of the |njf) states is numerically impossible as the cor- 
responding Hilbcrt space is of infinite dimension. How- 
ever, this problem can be circumvented by considering 
the rii^-resolved density matrices (jir S Z), which 

are submatrices of the complete density matrix p, whose 
entries are defined by 

(fc|p("-)|j) = (A:,n^|pb-,n«) , (19) 

where the relation between tir and n'^ is given by charge 
conservation. At the DJQP resonance, we have n'^ = 
UR - 2 for (fc,j) = (1,-1), n'p, ^ nn + 2 for (fc,j) = 
(— 1, 1) are ~ riR otherwise4^ Note that we did not 
write out the oscillator degree of freedom explicitly in this 
matrix. Calculating the time evolution of these matrix 
elements according to Eq. piT)) . one finds 

= [£ - 1,, J+p - J^p] (20) 

As expected, the tunneling leads to a coupling between 
density matrices of different tlr. It is produced by the 
current superoperators describing the quasiparticlc and 
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the Cooper pair tunneling, which are defined as 
T,pP = Tr\1){2\p\2){1\ , 



(21) 

^cpP = (-11 (-11 + Ii) (-11 HI) (1|] , 

|-i)(i|H-i)(-i| + |-i)(i|Hi)(i|] • 



'-cpP 



ijR 
2h 

ijR 
2h 



It is important to reahze that by writing the Liouville 
equation in terms of riij-resolved density matrices and 
current superoperators, we have achieved a description 
of the system in terms of only the |n,nosc) states again. 
This, however, comes at the price of having to deal with 
an infinite number of density matrices, p^'""-\ Still, fol- 
lowing the approach of Ref. [s^l it will turn out that 
convenient expressions for the current and the noise can 
be formulated in terms of these current superoperators. 



III. OSCILLATOR PROPERTIES 

As mentioned before we treat the NR harmonic 
oscillator and we use the master equation to investigate 
the time evolution of the mean displacement 



(x) = Tr„„ Tr„ Trosc [p(i)a 



(22) 



and of the velocity v, correspondingly. Here, Tr„(-) = 
J2n=-i ('^l ' I") denotes the trace over the island charge, 
while Tr„^ (-) = Y.nii=-oo ("-fli ' traces over the tun- 
neled charge. Likewise, the master equation will allow us 
to calculate expectation values of higher order like (x^) 
and (w^) which are required for the calculation of the 
oscillator energy. 

For a linear coupling of the NR to the SSET as in 
Eq. ([7]), we find the following equations describing the 
time evolution of the oscillator coupled to the SSET 

I (^) = («) , (23) 
v) = -n^ (x) - 7e.t {v) + 4 (") ' (24) 



dt 



M 



where (n) — k (pkk) is the expectation value of the is- 
land occupation and 7oxt accounts for the external damp- 
ing. The stationary Hmit, where -^(v) = -^{x) = 0, can 
be regarded as the long-time limit when the oscillatory 
behavior has been damped by the thermal bath and thus 
(v) = and (x) = {A/Mn'^){n). If the coupling A to 
the SSET is zero, the oscillator stays in its equilibrium 
position at (x) = 0. For finite coupling, due to the elec- 
tromagnetic repulsion the NR equilibrates in a position 
(x) ^ proportional to the coupling and the charge (71} 
on the island. Note that the influence of the SSET on the 
NR is of the first order in the coupling A. This regime has 
already been studied in some detail^*^ and it was shown 
that the SSET acts as an effective thermal bath for the 
NR. As we will illustrate further on, the signature of the 
NR in the transport properties of the SSET is of second 



order in the coupling and is clearly visible in the current 
and the noise properties of the SSET. 

To study the infiuence of the NR on the SSET we in- 
troduce dimensionless quantities which are normalized to 
motional quanta of the oscillator. Using the frequency f2 
and the harmonic oscillator length. 



2Am 



(25) 



as units we define x = x/xq, t = and v = v/VIxq^ 
i.e. we normalize all variables with respect to "oscillator" 
quantities. This allows an easier comparison with the 
experiment where, for example, the bias voltage can be 
varied at constant coupling. Consequently, the equations 
of motion can be rewritten as 



dt 
d 



dt 



-~{v) =- (i) - 7cxt (w) + 2A (n) 



(26) 
(27) 



where 7cxt = 7cxt/f^ and A = xqA/HQ. 

Not only is the equilibrium position of the resonator 
shifted by the coupling to the SSET, but also the cumu- 
lants of the position and velocity of the NR, i.e. ((x^)) = 
{x^) — (x)"^, are infiuenced by this coupling: 



f.{{{i,i}+))=2{{v'))-2{{£')) -7cxt(({i,«}+)) 
-h4i((ni) - (n) (i)) , 



(28) 



(29) 



dt 



r((«')) = - 27cxt((^')) +47cxtr£ 



+ AA {{n v) - (n) (v)) 



(30) 



where Tg = kBTg/hrt and {•,•}+ denotes the anti- 
commutator. In the stationary limit this leads to 

(({5,5}+)) -0, (31) 
{{v^)) - 2fB +2A{{nv)^ {n) {i)) /%^t , (32) 



= m+2A{{n S:)-{n) {£)) 



(33) 



To lowest (linear) order, which we refer to as the 
thermal-oscillator approximation, we assume that (nv) = 
(n) (v) . This is identical to assuming that the correla- 
tions between n and v vanish, i.e. {{nv)) = 0. Con- 
sequently, the fiuctuations of the harmonic oscillator 
are not influenced by the SSET such that the virial 
theorem ((w^)) = {{x'^)) and the equipartition theorem 
((w^)) ^ 2Tb are fulfilled in the high-temperature limit 
Tb 2> 1. The resonator is thus in a thermal state de- 
termined only by Tb and 7oxt- In the thermal-oscillator 
approximation analytic expressions for the current and 
noise in the SSET can be derived and will be discussed 
in the upcoming sections. 
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The thermal-oscillator approximation is justified for 
weak coupling between the SSET and the NR, but fails 
for stronger coupling, since the observablcs of the oscilla- 
tor become entangled with the charge state of the SSET. 
As was already observed before^^i^ an increased coupling 
can drive the oscillator to a non-thermal state character- 
ized by a finite (nv) — (n) (v) ^ 0, where the virial and 
equipartition theorems no longer hold. 

In order to investigate this regime, we have to go to 
the next order in our approximation which means taking 
the fluctuations of {{nx)) into account, but assuming all 
higher-order cumulants to vanish, e.g. {{nx^)) ~ 0. This 
will be referred to as the Gaussian approximation since 
for a Gaussian distribution all cumulants ((a;")) for n > 
2 are zero and the resonator is fully described by the 
two lowest moments. Under this assumption, we can 
express expectation values of the form {x^pkj ) as products 
of the lower-order expectation values {xpkj), (a;^), (x) 
and (pkj) ■ While this approach leads to a closed set of 
differential equations, the set will now be non-linear and 
has to be solved numerically. 

In principle, this approximation scheme can be contin- 
ued to even higher orders 4£ However, since the Gaussian 
approximation works well for the low-coupling regime we 
are interested in, we do not go beyond it. Ultimately, for 
even stronger coupling, the linear coupling between the 
SSET and the NR itself becomes questionable. 

In order to investigate the oscillator state in more de- 
tail, we study the energy E = ^MQ'^x'^ + ^Mv'^, in di- 
mensionless quantities 

Previous work^i^ has focused on the fluctuations of the 
number of charges on the island, n, which in linear re- 
sponse can be described by an effective damping and ef- 
fective temperature. We will discuss this approach in 
more detail in the context of the charge noise. For an 
identiflcation of the oscillator state, though, we choose a 
different route and investigate the energy of the NR. In 
the stationary limit, using Eqs. ([5^ and ([55]) . we flnd for 
the energy 



^ \ ^ ^ 4 2{{nv)) +%^,{{ni)) 1 2 



A flnite (x) ^ provides additional potential energy, but 
the contribution is small, since (x)^ = 4A^(n)^. There- 
fore, it is the correlations of the entangled SSET-NR sys- 
tem contained in the second term, which have the poten- 
tial to drive the system out of a thermal state. 

The results for a calculation of the energy in the Gaus- 
sian approximation for a typical, experimentally relevant 
set of parameters are shown in Figs. [3] and 0] where we 
display the oscillator energy as a function of gate voltage 
Vg and the bias voltage V measured away from the reso- 
nance position. As the DJQP cycle contains two Cooper 
pair tunneling events, there arc two resonance conditions 
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FIG. 3: Oscillator energy in units of HQ in the Gaussian ap- 
proximation as a function of the gate voltage eVa/ti^l and the 
bias voltage eV/hO, where (0,0) denotes the resonance. The 
parameters used are — Tr — 10, Jl ~ Jr = 2, 7oxt = 
lO""*, Tb = 2.5 and A = 0.02. In the red-detuned area 
{V < 0, —V < Vg < V), cooling below the bath tempera- 
ture is visible (blue region). Driving can be observed in the 
blue-detuned case {V > 0, -V < Va < V). The highest 
energies are obtained in the yellow region. 




-10 10 
eVc/m 

FIG. 4: Oscillator energy in the Gaussian approximation in 
units of hQ for the Gaussian approximation as a function 
of the gate voltage eVa /Kl and the bias voltage eV/ Kl for 
increased coupling A = 0.03. Cooling and driving effects are 
increased as compared to Fig. O In the black area two stable 
and one unstable solution are found, i.e. the system becomes 
bistable. The area grows for stronger coupling. 



which have to be met and which can be controlled by 
adjusting the bias and gate voltages. 

The physical picture can be explained most clearly if 
we assume Vg = 0, which corresponds to a vertical cut 
in Fig. [3l If the system is blue-detuned from a resonance 
{V > 0), the tunneling Cooper pairs transfer a part of 
their energy to the oscillator in order to be able to tunnel. 
This leads to driving of the oscillator. On the contrary, 
for a red-detuned resonance {V < 0), the Cooper pairs 
can absorb energy from the oscillator, leading to cooling. 
A similar result was already found using a linear-response 
approach in Ref. Q . 

In the regime where both resonances involved in the 
DJQP cycle are blue-detuned {V > 0, -I^ < Vg < V), 
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FIG. 5: Distribution P(nosc) of the oscillator phonon 
number for different values of the bias voltage eV/hQ = 
{— 1, 0, 1, . . . , 7} from left to right along the x-axis. The pa- 
rameters used for this plot are A = 0.1, Tl = Tfl = 12, 
Jl = Jr = 2.5, 7ext = 0.001 and Tb — 3. For negative V and 
small positive V we find an exponential decay corresponding 
to a thermal state. For larger V > Q the distribution develops 
a peak at n 7^ which indicates a driven state. 



we find a particularly strong driving of the oscillator. In 
the white regions of Fig. [H energies of the order IQ^Ml 
(depending on the system parameters) are reached even 
for rather small coupling of the order A k, 0.02. The nu- 
merical solution of the Liouville equation reveals more- 
over that the resulting oscillator state is highly non- 
thermal, i.e. the distribution function of oscillator states 
P{nosc) strongly deviates from a Boltzmann distribution. 

This is calculated in Fig. [5] using the numerical ap- 
proach for different values of eV/Ml. We find an ex- 
ponential decay for < corresponding to the high- 
temperature limit of the Boltzmann distribution and a 
trend towards a driven state for V > Q. 

In the regime where both resonances are red-detuned 
{V < Q, —V<Vg <V), we find a cooling of the oscil- 
lator to temperatures well below the bath temperature. 
This shows up in Figs. [3] and [4] as the little triangular- 
shaped regions below the center, where the oscillator en- 
ergy drops below the value corresponding to the bath 
temperature. 

Due to the non-linearity of the master equation, more 
than one physical solution may emerge and we find that 
this is indeed the case in the sector where the NR is 
strongly driven. An analogous effect was found previ- 
ously for the same system at the JQP cycled and for a 
more general class of systems i^'^1^^'^^1^^'^'^1^^ We find that 
generally, the response of the system close to a DJQP res- 
onance is much more pronounced than at the JQP in the 
sense that quantitatively similar effects may be observed 
at much smaller values of the coupling. This agrees with 
the prediction^ that the backaction effects at the DJQP 
exceed those of the JQP by a factor (P/J)*. Therefore 
the DJQP is favorable from the experimental point of 
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FIG. 6: Oscillator energy in arbitrary units as a function of 
bias voltage eV/Ml and coupling A calculated in the Gaussian 
approximation. For < 0, backaction leads to cooling of the 
oscillator (blue). On the contrary, strong driving (red/yellow) 
can be observed for V > Q. Above a critical coupling ^(V'), 
the system enters a bistable region (black). The parameters 
are the same as in Fig. [3] and eVc/Kl — 0. 



view since achieving a strong coupling is challenging. 

A plot of the location of the bistabilities found in the 
Gaussian approximation as a function of the bias voltage 

V and the coupling strength A is shown in Fig.[6l Again, 
cooling of the oscillator is seen in the blue regions for 

V < 0, whereas driving happens for > as is depicted 
by red regions. Towards stronger coupling, both effects 
increase in magnitude. Two stable solutions appear only 
for a blue-detuned SSET and the voltage range where 
such an effect is visible grows with increased coupling. 
When increasing the coupling for a given voltage V > 
(which corresponds to a vertical cut in Fig. [6]), the sys- 
tem will evolve from a thermal state via the bistable state 
to a single driven state. On the contrary, an increase in 
voltage (corresponding to a horizontal cut) carries the 
system from a thermal state to a driven state, then into 
the bistable region. Beyond the bistable region, the sys- 
tem will fall back to the thermal state. Note that the 
effect of driving is much stronger than the cooling of the 
NR (cf. Figs. El and HI). 

We have confirmed the existence of bistability using 
the numerical approach by explicitly calculating the com- 
plete probability distribution P{x) [results not shown ex- 
plicitly]. In a thermal state, this distribution shows a sin- 
gle peak at X = 0. In a driven state, on the contrary, two 
symmetric peaks at finite values \x\ ^ appear. In the 
bistable regime, the system switches between these two 
states, leading to a distribution P{x) that shows three 
peaks, one (thermal) at the origin and two side-peaks 
(driven) at \x\ > 0. From these studies of the bistable 
regime using the numerical approach, we noticed that 
the parameter range in which the system exhibits bista- 
bility is smaller than the one obtained via the mean-field 
solution. The structure of the bistable region could be 
better predicted using the analytical approach by includ- 
ing higher than second order cumulants, i.e. extending 
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the analysis beyond the Gaussian approximation. 

In the following sections, we will show that these dif- 
ferent states of the NR also manifest themselves in the 
transport properties of the SSET, i.e. the current and the 
current noise. 



IV. CURRENT PROPERTIES 

We showed in the previous section that the coupling 
of an SSET to an NR can drive the oscillator into a non- 
thermal state and effect in cooling, potentially even cool- 
ing down close to the ground statei^ii^ In the following, 
we will study if and how it is possible to measure signa- 
tures of the resonator state in the current and current 
noise characteristics of the SSET close to the DJQP res- 
onance. 

The number of electrons that have left the island to the 
right lead, n^j, is proportional to the transported charge 
and therefore determines the current flow. Hence, the 
expectation value of the current is given by 

{!) = {nn) = (-e) Tr{p(0 nn} . (36) 

Without loss of generality, we chose to measure the cur- 
rent across the right junction. In the stationary limit, the 
total current is conserved such that the currents across 
the left and right junctions are equal. In each DJQP 
cycle, two tunneling events take place at the right junc- 
tion, see Fig. [2] the transfer of a quasiparticle which 
takes the island from charge state |2) to the state 
Subsequently, a Cooper pair tunnels to the right lead 
and leaves the island in the state | — 1). Two processes 
involving only changes in ul and n, which therefore do 
not contribute to (/} , close the cycle in which 3 electrons 
in total have been transported through the island. 

Equivalently, in the stationary state the expectation 
value of the current (/) can be written using the super- 
operator formalism. From Eq. (|20[) . one finds 



{!) = (-e) Tr„ Trosc {ItotaiPstat) 



(37) 



where Itotai = ^qp — '^^cp '^-^cp ^^'^ superoperator 
describing the total current. Equation (|37p is used in this 
form in the numerical routine. 

For the analytic mean-field approximations we split the 
total current into two terms, (/) = (I^^) + (I^), cor- 
responding to a contribution from the tunneled Cooper 
pair (the non-dissipative current) , and a contribu- 
tion from the quasiparticle tunneling event (the dissipa- 
tive part) . For these two contributions we can write 
down the exact expressions for the dissipative 

(l^) = {-e)TT„ TrosciIqpPstat) , (38a) 
- {~e)TR {p22) , (38b) 
and non-dissipative part 

{I'''') = i-e) Tr„ Trosc {[2Icp - 2I+p]p,tat) , (39a) 

(39b) 



In the thermal-oscillator approximation, these expec- 
tation values can be calculated by solving for the cor- 
responding elements of the density matrix, as shown in 
detail in the Appendix |B1 We find that the vector of all 
finite ipkj) , |p), is given by |p) = iJiM.^^ \c) where M is 
the evolution matrix of the SSET system containing all 
the system parameters, Eq. (jB3[) . and the constant |c) is 
the inhomogcneous part of the master equation due to 
the normalization of the density matrix {pkk) = 1- 
Using this result the stationary current for the DJQP 
cycle can be written as 



(/) = -i~e)n 



1 



1 

r7 



1 



1 



(40) 



The inverse of the tunneling rates for quasiparticles and 
Cooper pairs, F^ jj and 7l._r, respectively, can be inter- 
preted as effective resistances for these processes. Then, 
Eq. (|40p is reminiscent of the current through a series 
of resistors, where the largest resistance determines the 
behavior. The Cooper pair tunneling rates are given by^ 



7L(i) = 2f 



(r^/2)2 + ei,o(i) 



7fl(i) = 2fL-~ 



•J s 



(Fi/2)2 + e2__^(j) 



(41) 
(42) 



where ekj denotes the difference in energy between the 
charge states \k) and \j) and thus measures the detun- 
ing from the DJQP resonance. The renormalized tun- 
neling rates of the SSET are defined by Fq. = Fa/fi and 
Ja ~ Ja/2Ml. If the Cooper pair tunneling, say, to the 
right lead is resonant, i.e. ei,-i = 0, the rate 7^ reaches 
a maximum at the value 7_r = SJ^/Fl = (2J^/Fi)/f2. 
It decays like a Lorentzian away from the resonance. Ex- 
pressions for the current in less general form are for exam- 
ple derived for e^-fc — in Ref. [43[ and for |ei__i| = |e2.o| 
in Ref. [s^. 

Due to the capacitive coupling of the SSET to the NR, 
the resonance is shifted compared to the uncoupled case 
in the thermal approximation. We find in dimensionless 
units (derivation given in Appendix |B]) 



eVr eV 
eVr eV 



(43) 
(44) 



(-e)2iJfl((pi,_i) 



where eV/hH and eVc/fi^ are the relative bias and gate 
voltages measured from the values at the DJQP reso- 
nance. 

The SSET is affected only if the average position of the 
NR is finite, i.e. (x) ^ 0. This shift in the equilibrium 
position of the NR effectively corresponds, from the point 
of view of the SSET, to a change in Vq and will therefore 
be referred to as an effective hackgate behavior later on. 
This effect is of second order in the coupling A since we 
observed in the previous section that {x) is linear in A. 
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FIG. 7: Difference of the inverse current (A)-l/(/> {A = 
0) versus bias voltage eVjM). for various coupling strengths 
A. From the intersection with the voltage axis, the shift 
2A (£) of the oscillator position can be read out. Lines cor- 
respond to the numerical analysis and agree well with the 
Gaussian approximation (points) in this parameter regime. 
Inset: Average current {/) through the SSET for the param- 
eters Vl=Vr = 10, Jl = Jr = 2, 7cxt = 10"", and Tb = 3. 



Note that the average displacement of the NR oscilla- 
tion in the stationary limit. 



(i) ^ 2A (n) ^ 2A 



1 



2(/) 
-1r{x) J 3(-e)ri 



(45) 



is determined by the Cooper pair tunneling rate 7^ in 
addition to the quasiparticle tunneling rate T p. This is in 
contrast to the JQP cycle where it is only the necessarily 
small l/Tfi which determines the displacement. Since 
the rates and the current are implicitly dependent on {x) 
via the Cooper pair tunneling rate, Eq. (|45p is a self- 
consistency equation. 

Going beyond the thermal-oscillator approximation, 
we use again the truncated master equation and the nu- 
merical approach to calculate the current via the general 
Eqs. ([55)) and ([M)) . In order to assess the quality of the 
Gaussian approximation, we first compare the results of 
the two approaches for low coupling strength. The re- 
sult is shown in Fig. [TJ where we plot the difference in 
the inverse current between the weakly coupled and the 
uncoupled system. The results of the Gaussian approxi- 
mation and the numerically evaluated lines arc in excel- 
lent agreement. The Lorentzian lincshape of the current 
(inset of Fig. ^ is preserved in case of the weak coupling. 

The change in the average current due to the coupling 
to the oscillator can be most transparently illustrated 
by plotting the difference of the inverse currents in the 
coupled and the uncoupled cases, see Fig. [T] As ob- 
vious from Eq. (j40p the inverse current 1 / (/) is given 
by the sum of rates involving the various transport pro- 
cesses. In the thermal oscillator approximation, the func- 
tion 1 / (J) {A) — 1/ {I) {A = 0) changes sign as a function 
of eV/Ml at a position which is proportional to 2A{x) as 




FIG. 8; Average current (7) versus bias voltage for A = 0.01 
and A = 0.032, Tb ~ 2.5, further parameters as in Fig.[7l For 
increased coupling the Lorentzian peak becomes distorted and 
two stable solutions emerge which result in two stable values 
for the current. 



is shown in Eq. (jBlSp . We expect this sign change to 
be the most feasible way to experimentally observe the 
influence of the NR on the SSET current and to investi- 
gate quantitatively the coupling strength using only the 
average current. 

The bistability of the oscillator states, which was al- 
ready discussed in the previous section, also manifests 
itself in the current through the SSET. For increased cou- 
pling, we find that the equation of motion derived within 
the Gaussian approximation has up to three solutions, 
of which two correspond to stable currents. The result- 
ing current-voltage characteristic is shown in Fig. [51 In 
the center, the usual DJQP resonance is clearly visible. 
While for very low coupling A w 0.01, the current still 
follows approximately a Lorentzian, strong deviations be- 
come visible already for A = 0.032. 

In an experimental setup, we do not expect two sta- 
ble currents to be distinguishable. Indeed, a current 
measurement of the SSET-resonator system will yield a 
weighted average^ because decoherence effects will lead 
to a switching between the two stable configurations on 
time scales large compared to the oscillation period but 
small compared to the measurement resolution)^ These 
switching rates can easily be inferred from a comparison 
of the measured current to the two stable values and from 
the current noise, as we shall show in section [Vl 



V. NOISE PROPERTIES 
A. Charge noise 

In the past, linear- response arguments^ have been used 
to support the idea that a generic detector acts on the 
resonator in the same way as a second thermal bath, and 
that the backaction on the resonator caused by charge 
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fluctuations on the island can be described essentially 
by two parameters, a damping rate 7eff and an effective 
temperature ksTee- These are related to the charge- 
fluctuation spectrum 



Sniuj) 



dt e'''^*((n(t)n(0))) 



/, N ~,2 Sni0j) - Sn{-Uj) 



27off('^) 



(46) 

(47) 
(48) 



in the limit ksTcffiuj) 3> huj^ where uj is given in units of 
left = kBTcs/hfl and 7cff = 7eff/fi. Since these ex- 
pressions follow from a linear-response calculation, both 
effective quantities arc written in terms of the bare 
charge-noise, calculated in the absence of coupling with 
the oscillator. 

Investigating the retarded and advanced (absorption 
and emission) contribution of the charge correlation ex- 
plicitly, we can derive an analytic expressions for <S'„(w) 
for the uncoupled SSET (see Appendix ID|) 



S„{uj) + Sni~Uj) 



(49) 



2(n|M[c.2 + M2] '(K+-(n))|p) , 

(50) 



where M denotes again the evolution matrix of the SSET 
system, {n\ is defined such that (n|p) = (n) and K± 
denote coupling matrices given explicitly in Eqs. (jCSp 
and jCEl). Note that ((n^)) = (n^) - (n}^ = (n| K+ |p) - 
(ji)^ and K_ acts only on the off-diagonal elements of 
(pkj) , i-fi- the Cooper pair tunneling terms. 

The self-consistency equation for ((i^)), that has to be 
solved in the Gaussian approximation, can be written as 



2^cxtTB + 27ssetJsset(£, ^^) 

7ext +7SSET(i,i^) 



(51) 



where 



7ssET(i,5') = i(2i)' (52) 
(n| (1 - 7ext (7oxt -f M)) (1 + 7extM + M2) "'k_ |p) , 

27ssETTssET(i, i') = {2A)^ (53) 
{n\ (7e.t + M) (1 + 7e,tM + M^y' (K+ - (n)) |p); 

We observe that the mean-field equation in second or- 
der provides the same physics as linear-response the- 
ory, i.e. 7cff at uj/fl = 1 is of the same form as 7sset- 
Since Eq. (jSTj) is a self-consistency equation for ((i^)) 
and not for the effective oscillator energy as in Eqs. ([47]) 
and (|48|) the expressions differ by a factor of 4. The 
result in Eq. ([^T]) is more accurate in the sense that 



the parameters of the damped oscillator are involved: 
1 + 7extM + M2 = (1 - (7ext/2)2) ^ (M + W2)^with 

a renormalized frequency of fir = ^ffi — (7oxt/2)^ and 
additional damping due to 7cxt/2. 

Note that Eq. ([5T|) is a self-consistency equation for 
(i^) since |p) = \p{x,x^)) and it has to be solved to- 
gether with (x) = 2A{n). Even if it is assumed that 
p) = |p(0,2Tb)), the expression contains a correction 
due to the finite quality factor of the NR. 

Whereas the approach describing the detector as an ef- 
fective bath proved very successful in providing a simple 
physical explanation of experiments,— some of its short- 
comings have started to be identified in recent theoret- 
ical works j^ii^ For example, it has very recently been 
proposed^ that the signature of the oscillator in the 
charge noise spectrum of a generic detector is not the 
one of a thermal oscillator. In the light of these findings, 
the calculation of the full frequency-dependent charge 
noise spectrum of the SSET near the DJQP in the pres- 
ence of an oscillator becomes relevant, even more so since 
the charge-noise spectrum is an experimentally accessible 
quantity. 

As shown in Appendix [d1 it is possible to use the 
master-equation approach to derive formal expressions 
for 5'„(w), at least for weak coupling in the thermal- 
oscillator approximation. However, it turns out that 
these expressions are difficult to evaluate explicitly for 
stronger couplings (in the Gaussian approximation). On 
the other hand, the fully numerical approach presented in 
scction lTl Cl ean easily be adapted to allow the calculation 
of finite-frequency correlation functions of system (as op- 
posed to bath) operators, using the quantum regression 
theoremj^^i^ In the following, we therefore discuss only 
the charge-noise spectrum obtained numerically. Note 
that we verified that our algorithm reproduces accurately 
the known charge and position fluctuation spectra in the 
uncoupled (A ~ 0) regime. 

Figure [9] shows the symmetrized (in frequency) charge- 
noise spectra, 5^y™(w) = [S'„(w)-HS'„(-tj)]/2, obtained 
for different values of the bias- voltage detuning from the 
DJQP resonance. The oscillator state, i.e. thermal or 
driven, can be determined from the Fock space probabil- 
ity distribution in Fig. O The inset shows that the sig- 
nature of the oscillator in S!^"'^™{lu) appears prominently 
around the natural frequency of the oscillator. Away 
from to ^ fl, the charge spectrum is only weakly affected 
by the oscillator, since the coupling of the island to the 
resonator changes the effective biasing conditions of the 
SSET. 

The main panel of Fig. [5] shows the evolution of the 
charge-noise spectra when the system is taken from the 
"cooling" region (V = —1) through the resonance point 
{V = 0), to the voltage regime where the state of the os- 
cillator becomes highly non-thermal. Unsurprisingly, the 
overall signal around ft increases dramatically when the 
oscillator enters the driven regime, reflecting the overall 
increase in the magnitude of Sx{uj fl) when the oscilla- 
tors energy is increased. Associated with this increased 
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FIG. 9: Frequency-dependent charge noise of the charge n on 
the SSET, symmetrized contribution S^^™™(tj), versus uj/Q. 
for different values of the bias voltage eV/h^ measured from 
the resonance. Parameters are identical to Fig. [5] While there 
is no structure in the case 1/ = 0, we can distinguish if the 
oscillator is driven or cooled by the symmetry of the peak 
aX u = ^. Inset: 5'^''™™ (a;) in a larger parameter regime; 
numerical calculation A = 0.1 compared to analytical result 
i = 0. 



magnitude is an overall reduction of the lincwidth, that 
is again explained rather straightforwardly via the de- 
creased total damping rate in this region (7cft < 0). 

The most interesting observation to be made about 
Fig. [5] is perhaps the striking similarity between the spec- 
tra at y < and V i). In both cases, we find not only a 
resonance at the rcnormalized frequency of the oscillator, 
but also a dip at its bare frequency, exactly as derived 
in Ref. [b^I for a generic detector. Note that the renor- 
mahzed frequency Vtr{V) = i/fi^ - [(^gxt -I- 7cff(T^))/2]2 
depends on the detuning V. For F < we find 7eff > 0, 
such that the renormalized frequency fJ,. is smaller than 
Q.. In the region y > and 7cft < 0, the situation 
is reversed and the resonance appears at frequencies fi,. 
higher than O, while the dip is pinned. No structure is ob- 
served exactly at ^ = which relates to 7eff(V = 0) = 
since absorption and emission of energy from the SSET 
to the oscillator at the DJQP resonance is equal. 

The presence of a dip is not compatible with a 
purely "thermal" state of the oscillator, even in cases 
where the Fock state probability distribution function 
P{nosc) decays exponentially like in the fully-thermal 
case [cf. Fig. [5]. Not only does this result confirm that the 
simple model used in Ref. j6^] also applies to the com- 
plex SSET-resonator system, it also demonstrates that 
the "Fano-like" lineshape, where both a resonance and 
a dip appear in the spectrum of the charge noise, char- 
acterizes nicely the charge noise spectrum on both the 
"driving" and "cooling" sides of the resonance. 



B. Current noise 

We argued previously that the non-thermal oscillator 
states (i.e., cooling and driving of the resonator) manifest 
themselves in the cumulants of coupled system. There- 
fore, it is reasonable to expect that the coupling with the 
NR will modify the current fluctuations of the SSET. 

In the following, we will focus mainly on the Fano fac- 
tor, i.e. the current noise at w = which is easily ac- 
cessible by standard noise measurement techniques. The 
frequency-dependent current noise is given by MacDon- 
ald's formula^ 



Sj{co)= / dte-*(({/(i),/(0)})) 



(54) 



From this formula, we can deduce the following zero- 
frequency limit,— 

d 



(55) 



Thus we have to determine the long-time limit of 
^((n^(t))). Note that we assume for the derivation sym- 
metric capacitances such that Sll{^ = 0) = Srr{uj = 0) 
due to charge conservation4i Since we use the sym- 
metrized current noise, the Fano factor is defined as 
F ^ Si{uj = 0)/(-e) (/) without a factor of 2. 

In order to use Eq. (|20p to calculate numerically the 
current noise, we follow closely the approach presented 
in Ref. [53| . The original approach which applied to in- 
coherent processes can be generalized to the coherent 
Josephson tunneling. The zero-frequency current noise 
is thus given by 



Siiiu ^ 0) = (-e) 



Tin Tfosc (-^no-isePsiai ) 
2 Tr^ ^"^oscP^total^^totalpstat) 



(56) 



where Ttotai — Iqp ~ ^" '^-^cp ^^"^ superoper- 



ator describing the total current as defined previously, 
liioise = Iqp + ~^ ^CP ^^'^ ^ ^^'^ pscudoinversc 

of the Liouvillian We compared both the Fano factor 
and the frequency-dependent current noise in the uncou- 
pled case with the exact expressions of Ref. IH (see also 
App.[E| and verified thus the correctness of Eq. ([55]) . 

As in the previous section it is possible to obtain an an- 
alytic expression in the thermal-oscillator approximation. 
For details of the calculation we refer to AppendixlE] We 
find 

^ _ ^ _|_ ^ ~/sym + /asym ^^^-^ 



1 



1 



1 



Tr IL IR 



The term 















\1L 





(58) 
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7 - FiA) - F(0) [A^] 
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FIG. 10: Difference in Fano factors AF = F{A) - F{A = 0) 
in units of as a function of bias voltage eV/Kl for Vl ~ 
Vr = 10, Jl = Jr = 2, 7ext = 0.0001 and Tb = 2.5. Lines 
depict results of the numerical approach. The results of the 
Gaussian approximation (points) fit quantitatively only for 
A < 0.01. Inset: the Fano factor F{A = 0.02) decreases at 
V f» when the Cooper pair tunneling is resonant. 

reduces the noise^i^ at the resonance; this is a conse- 
quence of the coherence of the Cooper pair tunneling. 
The asymmetric part 











- M 








\1L 


IRJ 



\lL 1r) {fRl2YJl (fi/2)2j2 

increases the noise level, for example if the coupling is 
asymmetric, ^ F/j, or the system is tuned away from 
the resonance, i.e. £2,0 7^ or ei__i 7^ 0. 

The coupling of the SSET and the NR enters via the 
{x) dependence of the Cooper pair tunneling rates 7l {x) 
and ^r{x) defined in Eqs. ((4T|) and ((42|) . In the thermal- 
oscillator approximation the effect is thus of second order 
in A as observed already for the current. In analogy to 
the previous section the coupling of the SSET to the NR 
leads in the thermal approximation only to a shift in the 
resonance position of the SSET. Since the resonance con- 
dition is important for coherent Cooper pair tunneling, 
the shift manifests itself in a higher noise signal around 
F = 0, leading to a peak in AF = F{A 7^ 0) - F(i = 0), 
see Fig. [TOl At large voltages both the Fano factor of the 
coupled and uncoupled SSET converge to the value 3/2. 

In order to access the regime of increased coupling, 
we investigated the current noise also in the Gaussian 
approximation and compared the result to the numerical 
calculation explained below. For very low coupling up 
io A = 0.01, as long as the resonator remains close to 
its thermal state, the two approximations coincide. For 
increased coupling, however, deviations from the effective 
backgate behavior appear. A plot of the difference in 
Fano factors AF as a function of the bias voltage is shown 
in Fig. \M 




eV/TiSl 



FIG. 11: Logarithmic Fano factor ln_F and oscillator energy 
E/hQ as a function of bias voltage eV/hfl. The parame- 
ters are eVa/hn = 2.5, Fl = Tr = 12, Jl ^ Jr = 2.5, 
7cxt = 0.001, Tb ~ 3, and A = 0.1. The Fano factors shows a 
first peak when the system enters the driven state. At the on- 
set of the bistability, the Fano factor increases drastically on 
the logarithmic scale, which we attribute to telegraph noise 
due to switching between the two stable configurations. In- 
set: frequency-dependent current noise in the bistable region 
{eV/hQ = 12.8) which decays on the scale of the switching 
rates. 



We find that the central peak in AF is accompanied 
by two negative side peaks where the current noise of the 
coupled SSET is lower than the uncoupled value. In the 
cooling and driving regimes the SSET absorbs and emits 
energy to the resonator in order to move closer to the 
resonance position. As pointed out before, the noise at 
the resonance is reduced due to the coherent Cooper pair 
tunneling and therefore the cooling and driving mecha- 
nisms are responsible for the negative side peaks. As in 
the energy of the resonator, Fig.[6l the driving is stronger 
than the cooling, which manifests itself in the heights of 
the negative side peaks. 

Even before the bistability arises, the coupling between 
the SSET and the NR leads to a modification of the cur- 
rent noise, from which one can deduce if (and to which 
extent) the oscillator is cooled or driven. 

In the bistable regime, the master equation approach 
leads to two solutions for the noise corresponding to the 
two stable solutions for the current. We argued previ- 
ously that in this regime, a thermal and a driven solu- 
tion coexist. Experimentally, however, the two solutions 
will not be completely stable and the system will switch 
between the two states on a time scale much slower than 
the other system time scales. As was shown generally 
for bistable systemsj^i^ and was also recovered for the 
case of the JQF;^ the experimentally measured noise will 
then be dominated by the switching between the two sta- 
ble states and will then essentially be telegraph noise. 
This type of low-frequency noise is most clearly visible 
in the Fano factor which we have evaluated using the 
numerical approach in Fig. 1111 
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The Fano factor first peaks roughly at the transition 
between the thermal state and the driven state. It then 
descends to a local minimum. This behavior is correlated 
to the stationary value of the average position of the os- 
cillator, (x) , which acts as an effective backgatc on the 
SSET as discussed in detail in connection with Fig. [TOl In 
the driven state the transfer of energy from the SSET to 
the NR is maximal and therefore we observe a minimum 
in the Fano factor. Beyond the pure driven state, the 
system reaches the bistable region which clearly shows 
up as a drastically increased Fano factor on the logarith- 
mic scale in Fig. [TTJ We attribute this zero-frequency 
noise feature to a slow switching between the two stable 
configurations of the system.— 

The inset of Fig. [Tl] shows the frequency-dependence 
of the current noise at the voltage where the Fano factor 
is maximal. As expected for a bistable system the value 
drops fast from the supcrpoissonian value on a scale given 
by the sum of the two switching rates^, which are very 
slow compared to the other scales of the system. In the 
frequency-dependent current noise [not explicitly shown 
here] we find signatures at the frequency of the oscilla- 
tor and at higher harmonics, which appear as a thermal 
resonance on top of the uncoupled current noise of the 
SSET. 



VI. CONCLUSIONS 

In conclusion, we have provided a comprehensive treat- 
ment of a superconducting single-electron transistor ca- 
pacitively coupled to a nanomechanical resonator. As- 
suming a linearized coupling for typical system parame- 
ters, we have found that signatures of the mutual inter- 
action manifest themselves as well in physical quantities 
related to the oscillator as in the current and noise prop- 
erties of the transistor. 

It is known that this setup with the SSET close to the 
resonance condition for the double Joscphson quasipar- 
ticle process is especially suitable to cool the oscillator 
or drive it. We confirm this behavior and explain how 
the different approximations fail or succeed in capturing 
these features. We use two complementary approaches, a 
purely numerical solution of the Liouville equation and a 
mean-field analysis, which lead to identical results in the 
accessible parameter region. In addition to driving effects 
of the oscillator, our solution predicts the emergence of 
a bistable regime for increased coupling. 

Although the consequences of cooling and driving the 
resonator have been studied in previous work, the ques- 
tions of if and how these non-thermal states can be mea- 
sured in the transport properties of the coupled SSET 
close to the DJQP resonance were not discussed in the 
literature. This obvious gap is addressed in our work. 
We find that the current through the SSET gives a mea- 
sure of the displacement of the oscillator since, at the 
lowest-order, the coupling with the NR modifies the gat- 
ing condition of the SSET. 



Investigating the fluctuations close to the resonance 
frequency of the mechanical oscillator, we find pro- 
nounced effects. For instance, the frequency-resolved 
charge noise shows a sharp resonance/dip structure close 
to the resonator frequency, which allows to estimate the 
sign of the effective damping. Moreover, we have stud- 
ied the zero-frequency current noise, i.e. the Fano factor. 
Depending on the applied voltages, the SSET coupled 
to an NR results in an increase or decrease of the Fano 
factor compared to the uncoupled case. We could show 
that a decrease is related to driving or cooling of the 
NR. Furthermore we find that the Fano factor increases 
notably as soon as the system enters the bistable state. 
This feature can be used experimentally to pinpoint the 
bistable region. The switching rates between the two 
stable regimes can be identified from the frequency de- 
pendence of the current noise. 

In comparison to an SSET driven at a JQP resonance, 
the DJQP shows stronger effects for lower couplings. Ex- 
tending the parameter regime to strong coupling, differ- 
ent methods have to be applied and also the assumption 
of a linear coupling is challenged. This will be the subject 
of future work. 
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APPENDIX A: COUPLING OF SSET AND NR 

The coupling between the SSET and the NR originates 
from a capacitive coupling between the island with the 
gated resonator. The displacement x, i.e. the deviation 
from the equilibrium distance d, determines the capaci- 
tance Cn{x) 



Cn{x) = eo — — = C% — C% ^ 



X + d 



1 + x/d' 



(Al) 



and thus the charging energy Ec and the charge on the 
island uq 



Ecix) 



1 



2 Cl+Cr + Cg + Cn{x) 
1 



1 - 



C'l, x/d 
Cs 1+x/d 



(A2) 



no{x) ^ -{ClVl + CrVr + CgVg + Cn{x)Vn) 



C%Vn x/d 

— Tin 



e 1 -|- x/d 



(A3) 
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The linear expansion of the charging term of the Hamil- 
tonian Hq in x yields 



Ec{x)[n + noix)] ^E%{n + n'',;) 

+ 2E%{ 



(A4) 



"oj — — 



We neglect a shift of the charging energy of the order 
of EcC^ / Cj: since in typical experimental setups^^ this 
prefactor is much smaller than the term in the second line 
of Eq. (jA4p . This term would be proportional to x{n + 
no)^ and would describe a slightly different coupling. 

Since ng constitutes a constant energy offset it is ne- 
glected in the calculation. The interaction constant in 
H^j = ~A n a; is thus given by 



(A5) 



and can be tuned by changing the gate voltage Vn ■ Note 
that we assume an experimental setup where the charge 
density on the dot can be changed by a (plunger) gate Vq 
and there is an additional gate voltage from the oscillator 
Vat. 

The dimensionless quantity A in our notation is thus 



A 

m 



Xo 



C% cVn xq 

Cs hn d 



(A6) 



With the parameters from Rcf. [IQ], for instance, we can 
get some very rough estimates for the renormalized quan- 



tities of the SSET. We find that the parameters are of the 
order of 



r„ = 2.5-17, 
J„ = 1 - 2.5 , 

A = 0.001 - 0.01 



(A7) 
(A8) 
(A9) 



APPENDIX B: MEAN-FIELD EQUATIONS FOR 
THERMAL-OSCILLATOR APPROXIMATION 



Using the thermal-oscillator approximation in Eq. ([Tt 
for the quantities {pkj) we get a closed set of equations 
for 



|P) = ((P-i,-i) 

(P2,2) (P2,0) (P0,2) 



(Bl) 



where (poo) = 1 - - {pi,i) - {P2.2) is determined 

by the normalization of the density, Tr[p] = 1. |p) fulfills 
a matrix equation 



dt 



j|p) = -M|p)+zJi|c) , 



(B2) 



where |c) = (0, 0, 0, 0, 0, 1, —1)"^ and the evolution matrix 
M is defined as 



M 



/ . -ijfl ijR 

-iJr iJr Tl/2 + iei^^i 

iJr -ijR . Tl/I - iei,_i 

\-iJL -iJh 



-r, 



2iJi ffly2-t-ie2o 
-2«Jl . f h/2 - ie2o/ 



(B3) 



where . stands for the entry 0. As discussed in the main 
text we normalize the SSET quantities with respect to 
resonator properties, hence the quasiparticle tunneling 
Tq = Ta/^^j the Cooper pair tunneling Jq, = Ja/2hfl 
and the resonance energies 



El, 



£2,0 



4Bc 



no 



eV 



2A (i) 



K + i)-^-2i< 



(B4) 
(B5) 



are shifted due to the coupling to the resonator in the 
thermal approximation. As can be seen in these expres- 
sions, the resonator is only coupled to the off-diagonal 



terms in the density matrix p, which correspond to 
Cooper pair tunneling events. 

Some assumptions are made to simplify to expressions 
for €2,0 and ei,_i. We assume that the SSET is tuned to 
the DJQP resonance with nj^'"^^ = — 1/2 and therefore 
we find deviations from it by 



AEc^ 



DJQP 



rDJQP 



where Vq is the gate voltage measured from Vq 
(symmetric bias voltage assumed). Furthermore there is 
a fixed value of the bias voltage V^'"^^ which is chosen 
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such that 



^yDJQP 4^ 



(B7) 



and we denote any shift away from V"^^^ as Ay. Con- 
sequently the energy detuning is given by the gate and 
bias voltage as 



£2,0 



(B8) 
(B9) 



For convenience only, it is assumed in the calculations 
that AeVo/fi^ and AeV/Hi^ have the same prefactors 
(see Eqs. (|43)) and (|44)) ') and for abbreviation also A is 
mostly skipped. 

In the stationary limit, i.e. d/dt |p) = 0, we find im- 
mediately 



|p) = iJlM-^ \c) 
or written down explicitly 



(BIO) 



IP) 4 



1 



ijR 



1 

1 



1 



1 



2 



--|-iei,-l 
1 



iJl 



I 

2 1 



_ ■+iC2,0 

2 1 



(Bll) 



This result can be inserted in Eqs. ((38|) and ([39|) and we 
find the result for the current (/) , Eq. (PD|) . 

Note that {I^^) = 2(1^) follows immediately from 
the equation (jB2p 



- (Pii) = - Jfl - (j5_i,i)) + Tr (p2,2) (B12) 

in the stationary limit. This expression is due to the 
nature of the DJQP process and does not change if we 
go to more sophisticated approximations. 

The shift of the oscillator resonance is proportional to 
the charge on the island, 



(n) = 2 (|52,2) + -{p-i,-i) 
1 I \ , f I 1 

,7i? TrJ \"fR IL 



1 

7-R 



1 



\1R 
2 {I) 



1 



1 



(B13) 



As obvious from Eq. (|40p the slowest rate limits the cur- 
rent and thus determines its value. In the symmetric 
case Jl ~ Jr = J and Fl = Vr = F the current can 
also be written more explicitly in terms of the system 
parameters^ 



6FJ2 



2e2o + 2e2_^-)-8J2-HF2 



(B14) 



As discussed in the main text the shift due to the coupling 
to the leads can be read off from 



3(-e)f7/ 1 1 

2 yji) " (/) (i = 0) 
1 1 



(B15) 



1 



1 



1l{x) "fLix = 0) Jr{x) 7_R(i = 0) " 

In the symmetric case Jl = Jr. = J and F^ = F^^ = F 
as assumed in the main text, this can be written more 
explicitly as 



3(-e)r2/ 1 1 

2 yji) " (/) (i = 0) 



2{2A{i))[2A{i)~2^-^ 



(B16) 



Since (x) depends on the bias voltage V ^ the value of 
2A{x) {V, Vg) can be extracted from the intersection with 
the voltage axis when plotting I /{I) [A) - l/(/)(A = 0). 



APPENDIX C: MEAN-FIELD EQUATIONS FOR 
GAUSSIAN APPROXIMATION 



In the thermal-oscillator approximation we assumed 
that ((nx)) = {{nvjj = 0. For stronger coupling, the ac- 
curacy of this assumption get worse and we proceed to 
calculate {{nx)) = (n||xp)) with (n| = (-1,1,0,0,2,0,0) 
in the Gaussian approximation. Here ||a;p)) denotes 
the vector of cumulants (^{xpkj)) = (xpkj) — (x) (pkj) ■ 
Higher-order expectation values are assumed to vanish, 
for example from {{x'^pkj)) = it follows that 

{x^Pk]) = 2 {x) {{xpkj) - {x) ipkj)) + (x^) {pkj) ■ 

(CI) 

This approximation leads us to the following set of equa- 
tions: 



-J |p) = -M |p) + zJl \c) + 2ziK_ I lip)) , (C2) 
dt 



d 



dt 

d_ 

dt 



. ||xp)) ^ -M ||ip)) + II^^P)) + 2iA((£^))K_ |p) , 



l-Dp)) = -M ||5p)) - ||ip)) - 7oxt ||wp)) 
-f2i(K+-(n))|p) , 



(C3) 



(C4) 
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where the couphng matrices are given by This set of equations can be solved in the stationary Hmit 

and is given here in matrix representation 

K+ =diag (-1,1, 0,0, 2,1,1) , (C5) 
K_ =diag (0,0, 1,-1,0, 1,-1) . (C6) 



|p) = Ul [m + i (2i)2K_ (1 + 7e,tM + M2) ' { (K+ - (n)) + i{{i')){%..t + M)K_ }] |c) 
pp)) = 2i(l + + m^y^ [(K+ - {n)) + i((5;-'))(7ext + M)K_] |p) . 

Since in the stationary hmit 

{n\\vT?)) = (n|M||Jp)) 
we do not need to calculate ((nS)) and find instead 

((i^)) = + 2i(7ipp)) = 2fB + 2i/7ext(n|M||ip)) + 2i(n||ip)) . 
Using the general solution for ||ip)) we can rewrite this equation with the expression 

270x1^5 + 27ssETJssET(i,i^) 



7oxt +7ssET(i,i^) 



where we define 



7SSET(i, i') = i{2Af{n\ (1 ~ 7ext (7cxt + M)) (l + 7extM + M^) 'k_ |p(i, i^)) 

(K+-(n))|p(i,i2)) . 



27ssETrssET(£, i') = (2A)2(n| (7^x1 + M) (l + 7extM + M2 

I 



(C7) 
(C8) 

(C9) 

(CIO) 

(Cll) 



The equation for ((i^)) has to be solved together with the 
expression for {x) which leads in some parameter regime 
to more than one solution. This bistability of the model 
is discussed in detail in the main text. 



APPENDIX D: CHARGE NOISE OF THE SSET 

To calculate the charge noise in the DJQP we evaluate 
the expression 



Sn{uj) 



dte"^\{n{t)n{0))) (Dl) 
dtcos{ujt){{{n{t)n{0))) + {{n{-t)n{0))) 
dts\n{Ljt){{{n{t)n{Q))) - ((n(-t)n(O))) 



+ i 



and take into account that Sn{t) is not necessarily sym- 
metric. Therefore we introduce two new functions^ 



X{t) = Tro3cTr„,,{e-^^*/''npe'«*/n , 



r]{t) = TroscTr„^{: 



,-iHt/h 



pne 



Ht/h 



} 



(D2) 
(D3) 



The partial traces are taken over the oscillator degree of 
freedom and the tunneled charge npi. These functions 



have the properties ((n(t)n(O))) ~ Tr„{nx(t)} — {n)'^ and 
((n(-t)n(O))) = Tr„{m7(t)} - (n)^ where Tr„ denotes the 
trace of the charge states of the island. The functions 
x{t) and ri{t) fulfill the same differential equation 



j.\x)^-m\x)+iJL{n) \c) 



dt ' " 



'M\t]) + Ul (n) |c) , 
with different initial conditions 



(D4) 
(D5) 

(D6) 
(D7) 

Straight- 

-(«))|P) 

((n(i)n(O))) - ((n(-t)n(O))) = 2 (n| e-^*K_|p) , (D8) 
and consequently 

M 



|X(0)) =Tro,cTr„,{np} = KJp) , 
|7?(0)) =TroscTr„, W = K,,|p) , 

where = K+ + K_ and = K+ - K_. 
forwardly we find 

((n(t)n(O))) + ((n(-t)n(O))) = 2 (n| e-^'"(K^ 



S'„(w) 



:2(n| 
2i (n\ 



M 



(K+-(n))|p) 
2K_|p). (D9) 



18 



For example the effective damping and temperature in 
linear-response theory is related to 



5„(w) + S'„(-tj) 



= 2i{n 



1 



+ M2 

M 



-M2 



K-|p), (DIO) 

(K+-(n))|p) , 
(Dll) 



as is discussed in more detail in the main text. Note the 
striking similarity to the self-consistency equation for the 
Gaussian approximation, Eqs. (|51II54|) . 



APPENDIX E: CALCULATION OF THE 
CURRENT NOISE 

We want to calculate d/dt {{n'j^{t)^) and study there- 
fore the time-evolution of 



d/dt (nlit)) ^Tr{p{t)<} 



(El) 



■ Tr { {priR + niip) [ur, H] ] 



Therefore we introduce the rtij-resolved density matrix. 
Defining the operator 



Pk 



\nR) {riR] 



(E2) 



such that {pkj) = Trosc Tr„^ Tr„{p|j)(fc| ® |nij;)(n'^|} and 
{riRPki) = TroscTr„„Tr„{pni{|j)(fc| ® \nR){n'j^\}. Note 

that {pkjnR) is a different quantity since ur and j5^j^'"" 
do not commute. 

We assume the stationary limit in the sense that {pkj ) 
does not depend on time. Then, the rii^-resolved den- 
sity operator can be calculated in the thermal-oscillator 
approximation where the differential equation reads 



dt 



~ IriRp) (t) = -M \nRp) + \A) + (ur) iJl \c) , (E3) 



1^) 



~2iJr \ 

-2ijR (pij) 
2ijR _i) 





V / 



(E4) 



since {—e){nR) = {I)t is explicitly time-dependent (in 
the stationary limit (/} {t) = {I)), the differential equa- 
tion has to be solved explicitly and we find 



\\nR,p)){i)^e-^^\\nRp(t. = Q))) 



(E5) 



+ (1 



-Mi 



1^) - 7^M-i |P) 
(-e) 



Note that In^p) denotes similarily to previous defini- 
tions the vector of finite {nRpkj) and \ \nRp)) ~ \nR,p) — 
{it-r) Ip). 

The noise itself has a structure similar to the current, 
but with some correction terms from the counting of ur 



(/) =(-e)(6|p) , (E6) 
d/di({nUi))) = 2{~ef{b\\nRp)) + {-efTR (p22> (E7) 

{b\ = (00 2ijR -2ijR f 0) . (E8) 



If we apply MacDonald's formula^ for the symmetrized 
noise, i.e. Si{t) = Si{—t), we find [neglecting a factor e^] 



Sdu^)-S,i.^0)^2it\M-^^-l^l^\A) 

-2(6|--fe#^Kp(t = 0)) , (E9) 



(tj/0)2 -hM2 

{I) 



Siiuj^O) =2{b\M-' \A) 



i-e) 



(ElO) 



+ f 7? {P22) - -ajR {{pi.-i) + (p-hi)) 



The explicit expression for Si{uj — 0) is given in the main 
text. 
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